Two qubits entanglement dynamics in a symmetry-broken environment 
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We study the temporal evolution of entanglement pertaining to two qubits interacting with a 
thermal bath. In particular we consider the simplest nontrivial spin bath models where symmetry 
breaking occurs and treat them by mean field approximation. We analytically find decoherence free 
entangled states as well as entangled states with an exponential decay of the quantum correlation 
at finite temperature. 
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I. INTRODUCTION 

<l> , 

[T i I Since 1935 |l|,|2j, entanglement has been recognized as one of the most puzzling features of Quantum Mechanics. 
However, it is nowadays a widespread opinion that it also represents a fundamental resource for many quantum 
information protocols. As such, entanglement deserves to be analyzed in all respects. A primary concern is its 
robustness against environmental effects, and a supplied literature exists aimed at preserving entanglement coherence 
0, S S 0. S 0- More recently, attention has been devoted to the problem of thermal entanglement 0] i- e - 
quantifying entanglement arising in spin chains at thermal equilibrium with a bath. In this approach environment 
determines the temperature T to allow for a thermal distribution of system energy levels, while the detailed interaction 
between system and environment is not an essential part of the matter. The same is true also for those works that 
focus on entanglement decoherence 0,0] (also known as disentanglement 0,0). In this context the study of 
entanglement time behaviour is carried on with a master equation formalism and markovian approximation |15| or, 
, more generally, with arguments provided by spin-boson models. 

In the present pa/per we are going to envisage a novel approach to the problem along the line introduced for the 
first time in Ref.[lj| (a similar outline but supported by numerical means is also present in 0]). There, the authors 
considered a one spin system interacting with a fermionic environment endowed with a structure capable of symmetry- 
breaking (18j |. It was shown by analytical methods that coherence time increases as magnetic order enlarges or, in 
other terms, as temperature decreases. Here we extend this argument to a two qubits system plunged in a fermionic 
environment described by Transverse Ising model (TIM) and Ising model (IM) [Tq . We shall examine the time 
evolution of concurrence of the bipartite system [Tfj| , and find environment-limited concurrences as well as unlimited 
ones according to environment ordering level. 

The paper is organized as follows: in section II we introduce the model by referring to p"fj| and we revise some 
results. In section III we extend the model to a bipartite systems, and we present the results of paradigmatic cases 
in Sec. IV. Finally, Sec. V is for conclusions. Explicit calculations are reported in Appendices A and B. 



II. THE MODEL 



We consider the general scenario of a system and a bath described by hamiltonians H s and Hb respectively, and 
interacting through the hamiltonian H s b- The total hamiltonian is then H = H s + H s b +Hb, and the initial density 
matrix is assumed to be factorized, i.e., p(0) = p s ® Pb- We are looking for the time evolution of the reduced system 
density matrix p s (t); in particular we are interested in its off diagonal elements, the so called "coherences". If H 
doesn't depend on time the total density matrix will evolve accordingly to 

p(t) = e-* Ht p(0y m . (1) 

We can then obtain the reduced density matrix by tracing out the bath degrees of freedom in Eq. Q 

p s (t) = tr BP (t) . (2) 

We now follow the line sketched in Ref.|l6| to introduce the model for a spin system interacting with a spin bath. 
First of all, we assume the bath density matrix having a thermal distribution, that is ps = (e~ H]B / T )/Z , with T the 
bath temperature multiplied by the Boltzmann constant, and Z = tr (e~^ B / T ) the partition function. Furthermore, 
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we ask the bath hamiltonian to be a "symmetry breakable" one, that is endowed with phase transition in the degrees 
of freedom that provide the coupling with the system. The simplest hamiltonian with these requirements is a long 
ranged Ising Model-like one (IM). We add to it a transverse field to include a more general case in the analysis, 
dealing eventually with a Transverse Ising Model bath hamiltonian (TIM). The differences between the two models 
are minimal as coherence and entanglement is concerning and, in any case, we will be able to find results for IM in 
the limit of no transverse field for TIM. These peculiar environment hamiltonians will be studied through mean field 
approximation 18]. 



A. TIM-environment 



Let us consider N + 1 spin-|, and let be the a component (a — x, y, z) of the jth spin (j = 0, 1, . . . , N). The 
label j = refers to the system operators while j = 1, . . . , N to the bath operators. Furthermore, = (SJ ± iSj)/2 
are the spin flip operators, and |0) and |1) are the lower and upper eigenstates of S z . The following hamiltonians 
define the energy of the system, of the TIM-bath and of the interaction between them: 

H s = - Mo S z , (3a) 
HsB = -^=S^S z k , (3b) 

H B = - w (3c) 

k i : k 

where fig is the coupling constant with an external magnetic field parallel to the z axis, Jo, J are exchange coupling 
constants and w is the strength of the transverse field; they are all non negative constants. The indices of the sums 
run from 1 to N. Eq.^) describes a material in which spins compete to align along the positive direction of x axis 
or along i axis following a ferromagnetic behaviour; of course in the latter case the absolute direction of alignment 
is not important since the hamiltonian is symmetric in z-operators. We can notice that energy exchanges between 
system and bath are not included in the interaction hamiltonian; this will generate a pure dephasing dynamics, in 
which energy will be conserved, and temporal evolution analytically solved. 

The main difficulty with Eqs.Q is represented by the nonlinear term in Hb- For this reason it is helpful to 
approximate it with a mean field bath hamiltonian, as explained in |l6j : 

H% f ^-wJ^SZ- ^Jm J2 S k+ m2jN ■ ( 4 ) 

k k 

In the above equation to is the order parameter of the phase transition. Its absolute value ranges from to ^ as long 
as temperature ranges from the critical value T c = i to : the greater \m\ the larger the magnetic order of the bath 
along z axis. In the following we are going to consider only positive values for m since results are sign- independent. 
Everything remains true with the substitution to — > —to. This is a consequence of Hb z-symmetry, that is not 
lost in . The order parameter m is implicitly defined by the following self-consistent equation for the quantity 
= ±V w 2 + 4m 2 J 2 (also O's sign, written here for sake of precision, is irrelevant, for the same reasons of to's): 

— = tanh — — . (5) 

It is worth noting that from Eq.(|5j) we have — * J for T — * 0; furthermore, from the definition of 0, we can see it 
tends to 2m J in the limit of no transverse field (w — > 0). 

Together with Eq.JSJ) we must consider the following condition on the transverse field to obtain an ordered phase 
with TIM: 

J < tanh (-) . (6) 

This condition is not satisfied in the range of temperatures above T c ; for this reason the whole formalism we are using 
is valid only in the broken phase. 

With the linearized mean field bath hamiltonian it is possible to evaluate the coherence of the system (see Appendix 
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Sa(t) 



tr B 
1 



Z 



tr B 



-H^ f /T\ AH 



= S (0)r TIM (t) 



where if m/ = ff. 



-ffg , and 



/ tmJJf)\ 

\qVn) 



8 . / tmJJo 

i — sin = 

J V BViV 



(7) 



(8) 



Equation tells us that the time evolution of the off diagonal term of the system density matrix, responsible for 
the coherence of the system, is enclosed in the time behaviour of the complex valued factor r TIM (t). In particular, in 
order to find system decoherence, we ask whether and when this factor's absolute value goes to zero. In the limit of 
large N we can approximate it as: 



~ exp 



J%m 2 t 2 



Z!-l 

e 2 



We can see from @ that the system coherence decays exponentially with time. The coherence time is: 



01 



j n m V J 2 - e 2 



(9) 



(10) 



and increases as temperature decreases; for T — it is r = oo, and the system remains coherent. This is quite a 
counter-intuitive effect since collective quan tum properties of materials endowed with phase transition disappear as 
ordering increases (see for instance Ref. [2(J])- The factor t 2 in the exponent denotes the intrinsically reversible nature 
of the pro cess, in contrast to irreversibility introduced by markovian approach, and is closely related to the "Zeno 
effect" [2lj. In particular the periodicity of r TIM (t) in Eq.(JHJ) leads to the so called "recoherences" on a Poincare time 
scale. Decoherence takes place in the limit of an environment with infinite degrees of freedom; besides, the same limit 
is necessary to support the mean field theory approach we adopted. Thus in this context the limit TV — > oo has a 
double function: to take into account the decoherence process and to give a meaning to the mean field approximation 
written above. 

We briefly notice here that the factor r TIM (t) in Eq.@ is exactly alike to r as (t) of Eq.(32) in [lfj. But as far as 
that paper is concerning we must point out some inaccuracies: the final result (32) is correct, but the intermediate 
steps to find it are not. In particular the general formula (11) applies only if the 3x3 matrices Aj, commute, and this 
is not true when you look at Eq.(29) of that article. For this reason the intermediate formula (30) is wrong and the 
oscillations showed in Fig.l are not present. 



B. Limit of no transverse field: IM-environment 



In the limit of w 



we obtain from J2J the IM-hamiltonians which lead to 

J§t 2 



\ r iM{t)\ ~ exp 



(11) 



We can notice the same behaviour as for TIM-bath, but slightly more transparent: the coherence time is t im = 

Jo y T^fn? an< ^ i ts UIIuts are t Im =Tc ' > = an< ^ t im =QS> ~ 00 • We- note that coherence explicit dependence on bath 
coupling constant J has disappeared in this case; only interaction coupling constant Jo enters coherence expression 
when the bath is an IM-one. Otherwise the J coupling is indirectly present in (|llf) because it has a role in determining 
the order parameter m by means of Eq.(|5jl- 



III. THE EXTENSION 



In this section we extend results obtained in the previous one by considering a two qubits system, and studying the 
time evolution of their entanglement. We assume that the system qubits, labeled by 01 and 02, interact between them 
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and with environment, that is symmetry-breakable and modeled by TIM hamiltonians generalizing those of Eqs.fpjl: 

(12a) 



H 



H s — £0^01 ^02' 
Jo 



Hp 



-w 



EQX ^_ \ y QZ QZ 

k ~ N ' J ' 



(12b) 
(12c) 



In above equations £o represents the coupling constant between the qubits. We have discarded both local interactions, 
like that between qubits and an external magnetic field, and local couplings with environment degrees of freedom, a 
situation resembling a "collective" system-environment pairing Q. 

As a measure of entanglement between two qubits we adopt the so called "concurrence" |19|. which ranges from 
for separable states to 1 for maximally entangled states. The concurrence is given by: 



C = max {Ai — A2 — A3 — A4, 0} . 



(13) 



where Ai, A2, A3 and A4 are the square roots of the eigenvalues, in decreasing order, of the matrix R = p a p s - Here p s 
is the density matrix of the 2 system qubits, and p s is the "time reversed" matrix given by 



p s = (o%! ® CT02) P*s (°01 ® °02) 



(14) 



where ct's are the usual Pauli matrices. The symbol p* s means complex conjugation of the matrix p s in the standard 
basis 1 00), |01), |10>, |11). 

We assume that the qubits are initially decoupled from the environment, and the bath having a thermal density 
matrix pb — (e~ HlB ^ T )/Z . Therefore, we can write the whole state as: 



with a generic system pure state: 

I*) = a |00) + (5 |01) + 7 |10) + 6 |11) , 



Pb 



M 2 + l/3| 2 + l7l 2 + H 2 = i 



(15) 



(16) 



The steps to find time evolution of Ea. i|15|) are similar to those leading to Eq.JJJ (see Appendix 0, but now operators 
are represented by 4 x 4 matrices, being our system composed by two qubits. After mean field approximation |@J for 
the bath hamiltonian and some elementary algebra we obtain the reduced density matrix as: 



/ H 2 a*(3A*e-^° a*jA*e-^° a*5B* 



Ps(t)=tr B ( P (t)) 



a/3*Ae^° \(3\ 2 /3* 7 f3*SA*e^° 

aj*Ae^° (3j* M 2 j*SA*e^° 

a5*B p6*Ae-* it t° ~f6* Ae~i^ \5\ 2 



(17) 



where the coefficients 



.4 
B 



cos 



tmJJi 



e 

+ 1— sin 



tmJJi 



N 



/2tmJJ \ .6 . f2tmJJ Q \ 



1 N 



(18a) 
(18b) 



characterize the time dependence of the concurrence. From the above expression of p s (t) we can find the matrix R(t) 
and its eigenvalues, and from them, as explained, the final concurrence of the system. The complete expression for 
R(t) and for coefficients of Eq. l|18[) is given in AppendixlBl In the following we are going to consider some paradigmatic 
cases for the initial state (11511. 



5 



IV. PARADIGMATIC CASES 



Case 1 



Let us set a = S = in Ea. (|16[l for the initial state of the system. We obtain 
reduces to: 



R(t) 



whose square rooted eigenvalues are: 



This leads to the following concurrence: 



/0 0\ 

2|/3| 2 | 7 | 2 2/3*|/3| 2 7 

2/3 7 *| 7 | 2 2|/?| 2 | 7 | 2 

\ o o o o/ 



Ai = 2|/3|| 7 | , 

A2 = A3 = A4 = , 



/3 1 01) + -7 1 10) and i? matrix 



(19) 



(20a) 
(20b) 



(21) 



The entanglement results time independent, so the state does not perceive the presence of the environment. The 
reason is that is an eigenstate of the interaction hamiltonian and so it represents a decoherence free entangled 
state jjj. Since w is not present in the concurrence written above we know that the expression for the concurrence 
would be exactly the same for an IM-environmcnt. 



2. Case 2 

Now we set f3 — 7 = in Ea. (|16fl and obtain the state |^) = a |00) + 5 |11). The R matrix becomes: 



R(t) 



(\a\ 2 \5\ 2 (l 






\B\ 



\ 2aS*\S\ 2 B 
with square rooted eigenvalues in decreasing order: 













2a* \a\ 2 SB* 






\a\ 2 \6\ 2 (l + \B\ 2 ) ) 



Ai = \a\\5\{\B\- 
A2 = \a\\6\(\\B\ 
A 3 = A 4 = 0. 

^From Eas. ((T%|l . for large N, we get: 

\B\ w exp 

Then, by using concurrence definition and Eas. l|23|l . we arrive at 



-UlmH 2 



1), 
-1|), 



J 2 
92 



C T1M =2H|5||5| = 2|a||5|ex P 



2,2 



e 2 



- 1 



(22) 



(23a) 
(23b) 
(23c) 



(24) 



(25) 



The time behaviour of the concurrence just obtained is shown in Fig^for different values of the ratio We notice 
that in this case the qubits perceive the presence of the thermal bath, which spoils entanglement between them; in 
fact the initial state is no longer an eigenstate of the interaction hamiltonian. Only for zero temperature the order 
parameter reaches its saturation value and the concurrence remains constant. The behaviour is very similar to that 
of one qubit system coherence described by Eq.Q, but entanglement decoherence is exactly twice faster than one 
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FIG. 1: Concurrence versus scaled time Jot. Curves from the left to the right are for ^- — {.75, .50, .35, .25}. The value of 
other parameter are w — 0.1, J — 2. 

qubit decoherence. This result agrees with what found in [12J. Furthermore, together with the previous case, it falls 
within the general limitations represented by the Universal Disentangling Machine |13| . 
In the limit ai->0we obtain the concurrence for an IM-bath: 



C IM =2|a||J|exp 



-2 jh 2 - - m 2 



(26) 



Analogously to what already noticed for the single qubit coherence, in this limit the factor J disappears from the 
explicit concurrence expression. The only exchange coupling constant that enters in the decoherence time for the 
concurrence is Jq. 



3. Case 3 



If we set a = j3 = we obtain a product state = 7 1 10) + 5 |11) = (7 |0) + S |1)) |1), which trivially gives: 

R(t) = (0) C = . (27) 
In this case TIM hamiltonians are not able to induce entanglement between system qubits. 



4- Case 4 



If we set a = /3 = 7 = <5=iwe obtain again a separable initial state, but different from the previous one: 
I*) = § (| 00) + |01) + 1 10) + |11)) = ^= (|0) + |1)) -±= (|0) + |1)). In this case the R matrix is not trivial: 



R ® = 16 



/ l + \B\ 2 -2\A\ 2 e~ u ^ 
2B - 2A 2 e~ u Z« 



&6 
\\A 

>,\A 



2-2\A\e lt t° 2-2\A\ e 
2-2\A\ 2 e l ^° 2-2\A\ « 



U io 2B* -2{A*y e' lt ^ 

2 At^ 



\ 



it£o 



£0 



-u 6o 

1 + \B\ 2 ~2\A\ 2 e-^° J 



where: 



U io = (2A*e-i lt t° - (A* + AB*) e** 



V io = (2Ae'^° — (A + A*B) e^*«° 



(28) 



(29a) 
(29b) 



The concurrence is evaluable esplicitly, but the expression is too much cumbersome therefore not reported here. We 
only show in Fig[21its behaviour. The concurrence starts from its null value and increases because of the interaction 
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25 50 75 




FIG. 2: Concurrence versus the scaled time Jot. The plot shows the limiting role of decoherence (dotted line) that falls down 
exponentially, on entanglement (continuous line). The value of parameters are w = 0.1, £ = 0.3, J = 2, ^ = .25. 

between system qubits. If there wasn't disentanglement it would reach its maximum and decrease again giving rise 
to oscillations of equal amplitude. Nevertheless, the presence of environment alters this temporal behaviour damping 
the oscillations. For suitable values of coupling constants it can even prevent qubits from entangling at all. The 
interesting question of the maximal entanglement generation under dephasing processes arises naturally in this case 

V. CONCLUSION 

We have studied time behaviour of entanglement between two qubits dipped in a large symmetry- breakable fermionic 
environment, below the critical temperature T c . In the frame of mean field theory analytical results are provided for 
concurrence of the bipartite system, with temperature as a parameter of the problem. The hamiltonians involved 
in the discussion are those typical of Transverse Ising Models (TIM), capable of magnetic ordering under suitable 
conditions. To assign them a physical meaning we notice that, upon addition of a transverse field in H s , our model 
resembles an array of Rydberg atoms interacting with a cavity mode of the radiation field [lfij . Nevertheless such an 
assumption for H s makes the problem unsolvable by analytical techniques, and requires numerical investigation that 
we plan to accomplish in a near future. Beside that an important improvement would be to overcome mean field 
approximations adopted in the text, by including the effect of fluctuations, or by applying the spin wave approach 
22] to the bath. 

What comes out from the paper is quite a counterintuitive conservation of entanglement in a bath with strong 
interactions: the bigger the coupling strength (or the lower the ratio the longer the time qubits remain entangled 
(Ea. (|25|l and Fig.QJ). In some cases entangled qubits don't perceive environment at all, and the system state is 
a decoherence free one (Ea. (|21|l ). Several connections with results from the field of entanglement decoherence are 
provided. We believe our analysis can be useful to complete knowledge about entanglement dynamical properties. 

Acknowledgements 
APPENDIX A 
1. Exponentiation of suitable matrices 

Let us define a 2 x 2 traceless matrix A as 

A = (aa x + ba z ) = { * \ V (Al) 
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with a, b real coefficients. The exponentiation of A gives: 

e A = (coshg) 



(cos q) I + i ( ) .4 



with q — \J (a 2 + b 2 ). Therefore: 

tr (e A ) = 2 (coshq) ;tr (e lA ) = 2 (cosq) 
Let us extend these arguments to three matrices 1, 1Z, and 1' of the same form of A: 



tr 



e a e n e a ' 



= tr 



(cos a;) 1 + i I I i. 



(coshy)/+ ^!^,7, 



(cos z) I + i 



sm 3 



J' 



(coshy) 



Sill \ / IjQjUli y 
2 (cos x) (cos z) + i (cos z) I I I I tr (11V) 



x J \ y 

, . / tanh«\ / sinz\ 
-i (cosx) tr (111 1 ) - 



sm.T \ / smz 



\ y J \ z J ' \ x J \ : 

where x, y, z are respectively related to the elements of 1, 1Z, 1' as q was related to A. 



tr (11') 



(A2) 



(A3) 



(A4) 



2. Coherence Expression for TIM 

As an example of calculation we report the steps that lead to Eq. 10 . All other calculations are easier than this one 
and can be performed following the same line. 
The time evolution of the total density matrix is: 



e -m 2 JN/T f f 
P (*) = 7? { exp < it 



-^=S Z + 2mJ j ,S7 - ,rS! 



exp \ (I /T) {wS% + 2m J SI) exp \ -it £ 



S z + 2m J SI + wS 



First, the partition function results: 



Z = e~ m2jN ' T tr^v 
By virtue of equation (|A3|) we find 



(l/T)"£(wS% + 2mJS* k ) 



_ e -rn\JN/T 



n*r[c 



(wSl+2 m JSl)/T 



Z = e 





' e " 




■|cosh 


') 




2T 





(A5) 



(A6) 



(A7) 



Notice that the constant e -™ 2 - JN / T j n the partition function simplifies with that present in Eq. l|A5l) . 

Let us now study the time evolution of the operator Sq = |0) (1| that represents the off diagonal part of the density 
matrix: 



So(t) 



2 cosh 



2T 



-, —N 



= 5 -(0) 



2 cosh I — 



-i -N 



Y[tr B {e a e K e a '} 



where: 



Jo 



1 = t 

' \2VN 
K = (wSl + 2mJS z k ) /T . 

Jo 



2m J SI + wSi 



1' = -t 



2VN 



2m J S z k + wS% 



(A8) 

(A9a) 
(A9b) 
(A9c) 
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In order to use Eg . (|A4|I we evaluate the following quantities: 

t 



8 /tanhy\ 2T 

2V Vn \n 



(AlOa) 
(AlOb) 
(AlOc) 



and 



^^/x t ( mJJn ^9 

,r(RI) = 2f br- e 

V 7 2 I AN 2 \N 



Then, substituting these into Eq. l|A4|) and performing the product we obtain: 



TTtrjeWn = 2* (cosh^-Y [cos f ^ U £ sin f 



JY 



(Alia) 
(Allb) 
(Allc) 

(A12) 



We can recognize in the second member of Eq. ljA12|) the constant r TIM (t) defined in Eq.(JSJ); the absolute value of it, in 
the limit of large N, gives the result of Eq.@. The other quantities of the article come out with similar calculations. 



APPENDIX B 



Complete R matrix for TIM 



Let's begin with time dependent density matrix expression for TIM hamiltonians Ijl2|l . After mean field approxi- 
mation we obtain: 



Pit) = ^ t 



I \ e -it(H s +H sB +H™f) psG -H™ f /Tjt{H s +H sB +H™ f ) 



= -g CX P \ H 



L k 

exp 



]T ( (S& + + 2m J )Si + Y, "St 



o' s e W \(l/T)J2(^ + 2mJSI) 



(Bl) 



Where we've have set p' s — e lt ^ aS o iS o 2 p s e~ lt ^ aS oi s <a . 

The constants present in Eas. (|18|) are found by complex conjugation of the following quantities, evaluated in a 
similar manner as the one seen in Appendix lAl 



]+wSZ]JwSZ+2rnJS* k )/T-itK^L+2mJ)SZ 



B* = 



- JJtTfl | e it [(-^+ 2m ^^+«'^] e (^+2mJS^)/T e -it[(^+2mJ)S^+»SJ]| 



(B2a) 
(B2b) 
(B2c) 



After calculations it's an easy task to verify that A* = D* , and for this reason the constant D doesn't appear in 
Eqs.CH. 
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R(t) = ( £ ^ ] (B3) 



The matrix for TIM is: 

— ( 

fi _ ( \a\ 2 \6\ 2 (l + \B\ 2 ^ -2a*(3-/8*\A\ 2 e^ 2a* (3 | 7 | 2 A*e~^° - \a\ 2 7 *c5 (A* + AB*) e i"& \ 

\ Q! / 3 *l'5| 2 (^ + ^*5)e^ t ?''-2|/3| 2 7^^e-^*?o -2a/3* 7 *c5 1 A| 2 e 4 *«° + 2 |/3| 2 | 7 | 2 / 

ff = / 2a* |/3| 2 7 A*e-^*«° - |a| 2 /3*(5 (A* +AB*)ei"& 2a* |a| 2 <5£?* — 2 (a*) 2 /3 7 (A*) 2 e~^ \ 

2 ~l^ -2a(/3*) 2 <5|A| 2 e 4 *« +2/3* |/3| 2 7 \a\ 2 (3* 5 {A* + AB*) - 2a* \(3\ 2 1 A* ) ( ' 

( a 7 * |<5| 2 (A + A*B*)e^° - 2(3 | 7 | 2 8*Ae~^ -2a ( 7 *) 2 S \A\ 2 + 2/3 7 * | 7 | : 



3 V 2a<5* |5| 2 B-2/3 7 («5*) 2 A 2 e-^» 2/3 1 7 | 2 5*Ae~i it ^° — a 7 * |<5| 2 (A + A*B) eA lt ^° J ( ' 

( -2a/3* 7 *(5 |A| 2 e^ + 2 |/3| 2 | 7 | 2 |a| 2 7 *<5 (A* + AB*) e^o - 2a* (3 | 7 | 2 A* e -^*?o \ 

^ 4 = y 2\p\ 2 jSAe-i u ^ -af3* \S\ 2 (A + A* B) \a\ 2 \S\ 2 (l + |B| 2 ) - 2a*(3 1 S* \ A\ 2 ) (B<) 

From it we have extracted all particular cases treated in the text. 
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